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\q" 1 Introduction 
o 

^ The roots of the theory of nonhnear dispersive waves date back to hydrodynamics of the 19th 
^ ■ century. It was observed, both experimentally and theoretically, that, under certain circum- 
Q ■ stances, the dissipative effects in nonlinear waves become less important then the dispersive 
| 7 | ■ ones. In this way, balance between nonlinearity and dispersion gives rise to formation of stable 
, patterns (solitons, cnoidal waves, etc.). Driven by applications in plasma physics, these phe- 
^ I nomena were widely studied, both analytically and numerically, starting from the middle of the 
■ 20th century. The main mathematical break-through in the theory of nonlinear evolutionary 
>■ PDEs was the discovery of the phenomenon of their integrability that became the starting 
^ ' point of the modern theory of integrable systems with Korteweg-de Vries equation being the 
O first instance in which integrability appeared. But most evolutionary PDEs are not integrable, 
^ of course. As a powerful tool for numerical simulations, the method of kinetic equation has 
\Q I been developed in 1960-th and applied to many different types of dispersive evolutionary PDEs. 
O ■ The wave kinetic equation is approximately equivalent to the initial nonlinear PDE: it is an 
' averaged equation imposed on a certain set of correlation functions and it is in fact one limiting 
case of the quantum Bose-Einstein equation while the Boltzman kinetic equation is its other 



Oh- 



^ limit. Some statistical assumptions have been used in order to obtain kinetic equations; the 
limit of their applicability then is a very complicated problem which should be solved separately 
for each specific equation. 



^ , The role of the nonlinear dispersive PDEs in the theoretical physics is so important that the 
^ " notion of dispersion is used for "physical" classification of the equations in partial variables. 

■ " " ' On the other hand, the only mentioning of the notion "dispersion relation" in mathematical 
literature we have found in the book of V.I. Arnold pTj who writes about important physical 
principles and concepts such as energy, variational principle, the Lagrangian theory, dispersion 
relations, the Hamiltonian formalism, etc. which gave a rise for the development of large areas in 
mathematics (theory of Fourier series and integrals, functional analysis, algebraic geometry and 
many others). But he also could not find place for it in the consequent mathematical presenta- 
tion of the theory of PDEs and the words "dispersion relation" appear only in the introduction. 

In our paper we present wave turbulence theory as a base of the "physical" classification of 
PDEs, trying to avoid as much as possible specific physical jargon and give a "pure" mathemati- 
cian a possibility to follow its general ideas and results. We show that the main mathematical 
object of the wave turbulence theory is an algebraic system of equations called resonant mani- 
folds. We also present here the model of the laminated wave turbulence that includes classical 
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statistical results on the turbulence as well as the results on the discrete wave systems. It is 
shown that discrete characteristics of the wave systems can be described in terms of integer 
points on the rational manifolds which is the main novelty of the theory of laminated turbu- 
lence. Some applications of this theory for explanation of important physical effects are given. 
A few open mathematical and numerical problems are formulated at the end. Our purpose is 
attract pure mathematicians to work on this subject. 

2 General Notions 

For the complicity of presentation we began this section with a very brief sketch of the tradi- 
tional mathematical approach to the classification of PDEs. 

2.1 Mathematical Classification 

Well-known mathematical classification of PDEs is based on the form of equations and can 
be briefly presented as follows. For a bivariate PDE of the second order 



and three types of PDEs arc defined: 

• b"^ < Aac, eUiptic PDE: ip^x + i^yy — 

• 6^ > 4ac, hyperbolic PDE: t/j^^ — '^yy — x'^x — 

• 6^ = 4ac, parabolic PDE: iIj^^ ~ '^xyipy — ■?/' = 

Each type of PDE demands then special type of initial/boundary conditions for the problem 
to be well-posed. "Bad" example of Tricomi equation yipxx + "4^^ = shows immediately 
incompleteness of this classification even for second order PDEs because a PDE can change its 
type depending, for instance, on the initial conditions. This classification can be generalized to 
PDEs of more variables but not to PDEs of higher order. 

2.2 Physical Classification 

Physical classification of PDEs is based on the form of solution and is almost not known 
to pure mathematicians. In this case, a PDE is regarded in the very general form, without 
any restrictions on the number of variables or the order of equation. On the other hand, the 
necessary preliminary step in this classification is the division of all the variables into two groups 
- time- and space-like variables. This division originated from the special relativity theory where 
time and three-dimensional space are treated together as a single four-dimensional Minkowski 
space. In Minkowski space a metrics allowing to compute an interval s along a curve between 
two events is defined analogously to distance in Euclidean space: 



yy ~ 



its characteristic equation is written as 



dx 
dy 




ds^ = dx'^ + dy^ + dz^ - c^dt^ 
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where c is speed of light, x, z and t denote respectively space and time variables. Notice that 
though in mathematical classification all variables are treated equally, obviously its results can 
be used in any applications only after similar division of variables have been done. 

Suppose now that linear PDE with constant coefficients has a wave-like solution 

t) = Aexp2[/cx — 0t] or t) = A sin(A;x — 0t) 

with amplitude A, wave-number k and wave frequency 0. Then the substitution of dt = 
—i0, dx = ik transforms LPDE into a polynomial on and k, for instance: 

i^t + ai^x + Pi^xxx = 0{k) = ak - (3k^ , 

i)tt + oi^<Pxxxx = =^ 0^{k) = a^k^, 
- « Vx. + = ^ 0\k) = + 
where a and (3 are constants. 



Definition Real-valued function = 0{k) : d^0/dA;^ ^ is called dispersion relation 
or dispersion function. A linear PDE with wave-like solutions are called evolutionary dispersive 
LPDE. A nonlinear PDE with dispersive linear part are called evolutionary dispersive NPDE. 

This way all PDEs are divided into two classes - dispersive and non-dispersive |2j. This 
classification is not complementary to a standard mathematical one. For instance, though 
hyperbolic PDEs normally do not have dispersive wave solutions, the hyperbolic equation iptt — 
Ci^'4'xx + /^^"^ = hcLS them. Given dispersion relation allows to re-construct corresponding 
linear PDE. All definitions above could be easily reformulated for a case of more space variables, 
namely Xi,X2, Linear part of the initial PDE takes then form 



^dt dxi dxn 
and correspondingly dispersion relation can be computed from 

P{-i0,iki, ...,ikn) = 

with the polynomial P. In this case we will have not a wave number k but a wave vector 
k = {ki, kn) and the condition of non-zero second derivative of the dispersion function takes 
a matrix form: 

1^1 ^0. 

dkidkj 



2.3 Perturbation technique 

Perturbation or asymptotic methods (see, for instance, [3]) are much in use in physics and are 
dealing with equations having some small parameter e > 0. To understand the results presen- 
tated in the next Section one needs to have some clear idea about the perturbation technique 
and this is the reason why we give here a simple algebraic example of its application. The main 
idea of a perturbation method is very straightforward - an unknown solution, depending on e, 
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is written out in a form of infinite series on different powers of e and coefficients in front of any 
power of e are computed consequently. Let us take an algebraic equation 

x'^ -{3-2e)x + 2 + 6 = (1) 

and try to find its asymptotic solutions. If e = we get 

- 3x + 2 = (2) 

with roots x = 1 and x = 2. Eq.([l]) is called perturbed and Eq.(l2]) - unperturbed. Natural 
suggestion is that the solutions of perturbed equation differ only a little bit from the solutions 
of unperturbed one. Let us look for solutions of Eq.([T]) in the form 

X = Xq + 6Xi + e^X2 + ... 

where Xq is a solution of Eq.®, i.e. Xq = 1 or Xq = 2. Substituting this infinite series into 
Eq.([T]), collecting all the terms with the same degree of e and consequent equaling to zero all 
coefficients in front of different powers of e leads to an algebraic system of equations 

'e^ :a;g-3a:o + 2 = 0, 

: 2x00:1 - 3xi - 2xq + 1 = 0, 
e : 2xqX2 + xi — 8x2 — 2xi = 0, 
I. ■ ■ ■ 

with solutions 

xo = 1, xi = — 1, X2 = 3, ... and x = 1 — e + 3e^ + 

Xq = 2, xi = 3, X2 = — 3,... and x = 2 + 3e — 3e^ + ... 
Notice that exact solutions of Eq. ([T]) are 

x = ^[3 + 2e± VTTSeT^] 

and the use of binomial representation for the expression under the square root 

-f— ) 

(1 + 8e + Ae^)^ = 1 + (8^ + Ae^) + ^^^(8e + 4^^)^ + ... = 1 + 4e - 6e^ + ... 
gives finally 

X = -{3 + 26 - 1 - Ae + Qe^ + ...) = l-e + 3e^ + ... 

and 

x = -{3 + 26 + 1 + As -66^ + ...) = 2 + 3e-3e^ + ... 

as before. This example was chosen because the exact solution in this case is known and can be 
compared to the asymptotic one. The same approach is used for partial differential equations, 
also in the cases when exact solutions are not known. The only difference would be more 
elaborated computations resulting in some system of ordinary differential equations instead of 
Eqs.(l3]) (see next Section). 
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3 Wave Turbulence Theory 



3.1 Wave Resonances 

Now we are going to introduce the notion of the wave resonance which is the mile-stone for the 
whole theory of evolutionary dispersive NPDE and therefore for the wave turbulence theory. 
Let us consider first a linear oscillator driven by a small force 



Here p is eigenfrequency of the system, Q is frequency of the driving force and e > is a 
small parameter. Deviation of this system from equilibrium is small (of order e), if there is no 
resonance between the frequency of the driving force ee**^* and an eigenfrequency of the system. 
If these frequencies coincide then the amplitude of oscillator grows hnearly with the time and 
this situation is called resonance in physics. Mathematically it means existence of unbounded 
solutions. 

Let us now regard a (weakly) nonlinear PDE of the form 



where L is an arbitrary linear dispersive operator and N is an arbitrary nonlinear operator. 
Any two solutions of L^ip) = can be written out as 



with constant amplitudes Ai,A2. Intuitively natural expectation is that solutions of weakly 
nonlinear PDE will have the same form as linear waves but perhaps with amplitudes depending 
on time. Taking into account that nonlinearity is small, each amplitude is regarded as a slow- 
varying function of time, that is Aj = Aj{t/e). Standard notation is Aj = Aj{T) where T = t/e 
is called slow time. Since wave energy is by definition proportional to amplitude's square A'j it 
means that in case of nonlinear PDE waves exchange their energy. This effect is also described 
as "waves are interacting with each other" or "there exists energy transfer through the wave 
spectrum" or similar. 

Unlike linear waves for which their linear combination was also solution of L(^) = 0, it is 
not the case for nonlinear waves. Indeed, substitution of two linear waves into the operator 

— * — * — * — * 

eN{il)) generates terms of the form expi[{ki + k2)x — [0{ki) + uj{k2)]t which play the role of 
a small driving force for the linear wave system similar to the case of linear oscillator above. 
This driving force gives a small effect on a wave system till resonance occurs, i.e. till the wave 
number and the wave frequency of the driving force does not coincide with some wave number 
and some frequency of cigenfunction: 



where notation 0j = u!{ki) is used. This system describes so-called resonance conditions or 
resonance manifold. 



Xtt + P^x = ee' 



(4) 



Aiexpi[kix — u!{ki)t] and A2expi[k2X — U!{k2)t] 




(5) 
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The perturbation technique described above produces the equations for the amphtudes of 
resonantly interacting waves Aj = Aj(T). Let us demonstrate it taking as example barotropic 
vorticity equation (BVE) on a sphere 

^^^+2^ + .J(^,A^) = (6) 



dt dX 



where 



. , ^'^^p 1 d'^ip ^dip . T/ 7\ 1 ,dadh dadb, 
^^ = 7^^ + — ^T^^^^^^^T^ J{a,b) = -7^777-7777^7-- 

0(p^ COS^ (p OA^ 0(j) COSCp OA OCf) OCf) OA 

The linear part of spherical BVE has wave solutions in the form 

AP„"(sin0) expz[mA + , 

where A is constant wave amplitude, uj = —2m/[n{n + 1)] and P^{x) is the associated 
Legendre function of degree n and order m. 

One of the reasons to choose this equation as an example is following. Till now a linear 
wave was supposed to have much more simple form, namely, Aexpi[kx — u!{k)t] without any 
additional factor of a functional form. For a physicist it is intuitively clear that if the factor is 
some oscillatory function of only space variables then we will still have a wave of a sort "but 
it would be difficult to include it in an overall definition. We seem to be left at present with 
the looser idea that whenever oscillations in space are coupled with oscillation in time through 
a dispersion relation, we expect the typical effects of dispersive waves" p]. By the way, most 
physically important dispersive equations have the waves of this form. 

Now let us keep in mind that a wave is something more complicated then just a sin but still 
smooth and periodic, and let us look where perturbation method will lead us. An approximate 
solution has a form 

= M\ <P, t, T) + e^lj^{\, 0, t, T) + £2^2(A, 0, t, T) + ... 

where T = t/e is the slow time and the zero approximation ipQ is given as a sum of three linear 
waves: 

3 

MK 0, t,T) = J2 MT)P^^^ exp(«^fe) (7) 
fc=i 

with notations P^'^-' = P^'' and 6k = rn^X — uJkt- Then 



dA^o/dt + 2d^o/dX = 0, 

dAi/ji/dt + 2d'^i/dX = - J(V^o, Ai/jo) - dAiPo/dT, (8) 
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and 



I 



dX 

fe=i 



OA/ 3 7 



d<t> # 



J(V^o, AV^o) 



3 , 

-i iVfcm,p(^)— p(^)[A, exp(ie,0 -^,*exp(-i^,)][>lfeexp(i^fc) +A*exp(-i^^ + 
3,k=i "^"P 

3 , 

+t Y NkmkP^'^—P^^^[Ajexp{iej)+A*exp{-iej)][Akexp{^^^^ 
j,k=i ^ 

^ = i:/"'='A'.l^exp(,«.) + ^exp(-».)l. 

fc=l 

leads to the condition of unbounded growth of the left hand in the form 



J(^o, A^o) 



dT 



with resonance conditions 9j + 6^ — Oi V j, A;,i = 1,2,3. Let us fix some specific resonance 
condition, say, Oi + 02 — O3, then 

/(V-o, AV'o) ^ -i(JVi - JV2)(m2P<">-^P'") - miPW-^pW ■ 

>liA2exp[«(^i + ^2)] - ^tA;exp[-^(^i + ^2)], 

where notation = means that only those terms are written out which can generate chosen 
resonance. Let us substitute these expressions into the coefficient by s^, i.e. into the equation 

^ + 2^ = -J(^o,A^o)-^, 
multiply both parts of it by 

P(3) sin 0(^3 exp(i^3) + A* exp(-i^3)] 
and integrate all over the sphere with i — > 00. As a result following equation can be obtained: 

N,^^2zZ{N2-N,)A^A2, 

where 

J-n/2 d(l> d(t^ d(t> 
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The same procedure obviously provides the analogous equations for A2 and A3 while fixing 
corresponding resonance conditions: 



dAi 



2lZ{N2 - iV3)A3A; 



N2 



dT 
dA2 



2iZ{N3 - Ni)AlA: 



dT 



(9) 



In general, the simplest system of equations on the amplitudes of three resonantly interacting 
waves is often regarded in the form 



and is refereed to as a 3-wave system (keeping in mind that analogous system has to be written 
out for A* ). Coefficients depend on the initial NPDE. Similar system of equations can be 
obtained for 4-wave interactions, with the products of three different amplitudes on the right 
hand, and so on. 

3.2 Zakharov-Kolmogorov energy spectra 

The idea that a dispersive wave system contains many resonances and wave interactions are 
stochastic led to the statistical theory of wave turbulence. This theory is well developed [Ij 
and widely used in oceanology and plasma physics describing a lot of turbulent transport 
phenomena. Avoiding the language of Hamiltonian systems, correlators of a wave field, etc., 
one can formulate its main results in the following way. Any nonlinearity in Eq.(jl]) can be 
written out as 



where 6 is Dirac delta-function and V(i2..i) is a vertex coefficient. This presentation, together 
with some additional statistical suggestions, is used then to construct a wave kinetic equation, 
with corresponding vertex coefficients and delta-functions in the under-integral expression, of 
the form 



for 3-waves interactions, and similar for i-waves interactions. One of the most important 
discoveries in the statistical wave turbulence theory are stationary exact solutions of the kinetic 
equations first found in [5]. These solutions are now called Zakharov-Kolmogorov (ZK) energy 
spectra and they describe energy cascade in the wave field. In other words, energy of the wave 
with wave vector k is proportional to fc" with a < and magnitude of a depends on the specific 
of the wave system. Discovery of ZK spectra played tremendous role in the wave turbulence 
theory and till the works of last decade [Tj it was not realized that some turbulent effects are 
not due to the statistical properties of a wave field and are not described by kinetic equations 
or ZK energy spectra. 




(10) 




(11) 
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3.3 Small Divisors Problem 



In order to use presentation ( ITTl) one has to check whether so defined nonhnearity is finite. 
This problem is known as the small divisors problem and its solution depends on whether wave 
vectors have real or integer coordinates. 

Wave systems with continuous spectra were studied by Kolmogorov, Arnold and Moser [6] 
and main results of KAM-theory can be briefiy formulated as follows. If dispersion function 
is defined on real- valued wave vectors and the ratio a^j = 0i/0j is not a rational number for 
any two wave vectors ki and kj, then 

• (IC) Wave system is decomposed into disjoint invariant sets (KAM tori) carrying quasi- 
periodic motions; 

• (2C) If the size of the wave system tends to infinity, (IC) does not contradict ergodicity, 
random phase approximation can be assumed, kinetic equations and ZK energy spectra 
describe the wave system properly; 

• (3C) Union of invariant tori has positive Liouville measure and Q has measure 0, there 
exclusion of the waves with rational ratio of their dispersions is supposed to be not very 
important. 

Wave systems with discrete spectra demonstrate [7J substantially different energetic behav- 
ior: 

• (ID) Wave system is decomposed into disjoint discrete classes carrying periodic motions 
or empty; for any two waves with wave vectors ki and kj belonging to the same class, the 
ratio aij = 0i/0j is a rational number, 

• (2D) Energetic behavior of the wave system does not depend on its size, is not stochastic 
and is described by a few isolated periodic processes governed by Sys.lITOl): 

• (3D) In many wave systems (for instance, in the systems with periodic or zero boundary 
conditions) KAM-tori do not exist and discrete classes play the major role in the energy 
transfer. 

3.4 Laminated Wave Turbulence 

The results formulated in the previous section gave rise to the model of laminated wave turbu- 
lence [9] which includes two co-existing layers of turbulence in a wave system - continuous and 
discrete layers, each demonstrating specific energetic behavior. In other words, KAM-theory 
describes the wave systems leaving some "holes" in the wave spectra which are "full- filled" in 
the theory of laminated turbulence. 

Continuous layer, with its kinetic equation, energy cascades, ZK spectra, etc. is well-studied 
while the existence of the discrete layer was realized quite recently. In order to understand which 
manifestations of the discrete layer are to be expected in numerical or laboratory experiments, 
let us regard an example with dispersion function = 1/ ym^ + n? . First of all, it is important 
to realize: the fact that the ratio aij = 0i/0j is a rational number does not imply that dispersion 
function is a rational function. Indeed, for = 1/ \^m? + n^ we have 

: Z X Z ^ M 
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and for wave vectors fci = (2, 1) and ^2 = (9, 18), the ratio 0i/02 = 1/3 is rational number 
though is irrational function of integer variables. Decomposition of all discrete waves into 
disjoint discrete classes Clq in this case has the form 

{ki = {mi,ni)} e Clg if \ki\ = -fi^/q, W i = 1,2, ... 

where 7^ is some integer and q is the same square-free integer for all wave vectors of the class 
Clq. The equation 0i + 02 = 03 has solutions only if 

3g : ki, k2, k^ G Clg. 

This is necessary condition, not sufficient. The use of this necessary condition allows to cut 
back substantially computation time needed to find solution of irrational equations in integers. 
Namely, one has to construct classes first and afterwards look for the solutions among the waves 
belonging to the same class. In this way, instead of solving the irrational equation on 6 variables 

l/-^mf + n2 + l/^m2 + ni = l/^Jml + nl (12) 
it is enough to solve the rational equation on 3 variables 

1/71 + 1/72 = 1/73- (13) 

Some classes can be empty, for instance CI2 and C/3 in our example. Obviously, each non- 
empty class has infinite number of elements due to the existence of proportional vectors so 
that ki = (2, 1) with the norm = -\/5 and its proportional k2 = (9, 18) with the norm 
1^2! = 9v^ belong to the same class C/5. On the other hand, not all the elements of a class 
are parts of some solution which means that not all waves take part in resonant interactions. 

Let us come back to physical interpretation of these results. Resonantly interacting waves 
will change their amplitudes according to Sys. lITOl) . In this case the role of ZK spectra fc", a < 
0, is played by the interaction coefficient Z ~ A;", a > (see Fig.l). Non-interacting waves 
will have constant amplitudes (they are not shown in Fig.l). In the next Section we demonstrate 
some examples of different wave systems whose behavior is explained by the theory of laminated 
turbulence. 

4 Examples 

• Ex.1 Turbulence of capillary waves (dispersion function uj'^ = k^, three- wave interactions) 
was studied in [|10] in the frame of simplified dynamical equations for the potential flow 
of an ideal incompressible fluid. Coexistence of ZK energy spectra and a set of discrete 
waves with constant amplitudes was clearly demonstrated. The reason why in this case 
the amplitudes are constant is following: equation 

,3/2 ,3/2 _ ,3/2 

has no integer solutions [7] . It means that there exist no three- wave resonant interactions 
among discrete capillary waves, they take no part in the energy transfer through the wave 
spectrum and just keep their energy at the low enough level of nonlinearity. 
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Figure 1: Two layers of turbulence are shown symbolically. Low panel: 2D-domain in spectral 
space, nodes of the integer lattice are connected with the lines which correspond to 3-wave 
resonant interactions of the discrete layer. Middle panel: ZK energy spectrum with "the 

holes" in the nodes of the integer lattice. Upper panel: Interaction coefficient Z ~ /c^/^ in 
the nodes of the integer lattice. 



• Ex.2 Similar numerical simulations [TT] with gravity waves on the surface of deep ideal 
incompressible fluid (dispersion function uj'^ = k^, four-wave interactions) show again 
coexistence of ZK energy spectra and a set of discrete waves. But in this case waves 
amplitudes are not constant any more, discrete waves do exchange their energy and in 
fact play major role in the energy transfer due to the fact that equation 

1.1/2 ,1/2 _ ,1/2 ,1/2 

has many non-trivial integer solutions pKj (it is important in this case that 2-dimensional 
waves are regraded, i.e. k = \k\ = yrn? + ri? with integer m,n). 

• Ex.3 Some recurrent patterns were found in different atmospheric data sets (rawind- 
sonde time series of zonal wind, atmospheric angular momentum, atmospheric pressure, 
etc.) These large-scale quasi-periodic patterns appear repeatedly at fixed geographic lo- 
cations, have periods 10-100 days and are called intra-seasonal oscillations in the Earth 
atmosphere. In [13] Eq.([n]) (dispersion function u = —2m/n{n + 1), three- wave interac- 
tions) is studied which is classically regarded as a basic model of climate variability in 
the Earth atmosphere. It is shown that a possible explanation of the intra-seasonal oscil- 
lations can be done in terms of a few specific, resonantly interacting triads of planetary 
waves, isolated from the system of all the rest planetary waves. 

Remark In contrast to the first two examples, in this case only discrete layer of 
turbulence exists. Indeed, while = —2m/n{n + 1) is a rational function, any ratio 
0{mi,ni)/0{mj,nj) is a rational number and KAM-theory is not applicable. 

• Ex.4 A very challenging idea indeed is to use the theory of laminated turbulence to 
explain so-called anomalous energy transport in tokamaks. Turbulent processes respon- 
sible for these effects are usually described as H- and L-modes and ELMs (high, low and 
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edge localized modes consequently). Interpretation of the known experimental results in 
terms of non- resonant (H), resonant (L) and resonant with small non-zero resonance width 
(ELM) modes gives immediately a lot of interesting results. In this case non-resonant 
discrete waves are of major interest because they will keep their energy as in Ex.1, for a 
substantial period of time. This approach allows to get two kind of results: 1) to describe 
the set of the boundary conditions providing no resonances at all - say, if the ratio of 
the sides in the rectangular domain is 2/7, no exact resonances exist; or 2) to compute 
explicitly all the characteristics of the non-resonant waves (wave numbers, frequencies, 
etc.) for given boundary conditions. Some preliminary results are presented in [Hj, in 
the frame of Hasegawa-Mima equation in a plane rectangular domain with zero boundary 
conditions (dispersion function cu = \j^n^ -|- m^, three- wave interactions). 

It is important to remember that a choice of initial and/or boundary conditions for a specific 
PDE might lead to a substantially different form of dispersion function and consequently to 
the qualitatively different behavior of the wave system. For instance, Ex.3 and Ex.4 are 
described by the same Eq.(l6]) regarded on a sphere [rational dispersion function, only discrete 
layer of turbulence exists) and in a rectangular [irrational dispersion function, both layers exist) 
respectively. For some equations, a special choice of boundary conditions leads to transcendental 
dispersion functions. 

5 Open Questions 

We have seen that the main algebraic object of the wave turbulence theory is the equation 

0(mi, ui) + 0(m2, n2) + ... + 0(ms, Ug) = (14) 

where dispersion function is a solution of a dispersive evolutionary LPDE with | g^.^ | 7^ 0. 
So defined class of dispersion functions includes rational, irrational or transcendental function, 
for instance 

0[k) = ak - jSk^ , 0^[k) = a^k"^ + f3'^, = m/[k + l), = tanh afc, • • • 

where a and P are constants and k = y/m"^ + n^. Continuous layer of the wave turbulence, 
that is, with mj, Ui G K, is well studied. On the contrary, there are still a lot of unanswered 
questions concerning the discrete layer of turbulence, mj, rij G Z, and we formulated here just 
a few of them. 

Eq.f fT^ can be regarded as a summation rule for the rational points of the manifold given 
by 0. These manifolds have very special structure - namely, they can be transformed into a 
one-parametric family of simpler manifolds, namely fll2p into (|T^ . This situation is general 
enough, the definition of classes can be generalized for a given c G N, c 7^ 0, 1, — 1 considering 
algebraic numbers /c^^'^, /c G N and their unique representation 

kc = 7 e Z 

where g is a product 

while pi, ...pn are all different primes and the powers ei, ...e„ G N are all smaller than c. Then 
algebraic numbers with the same q form the class Clq and the following statement holds: 
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The equation aiki + a2k2... + a„A;„ = 0, Oj G Z where each ki = 'yiq^^ belongs to some class 
Qi ^ 1i)Q.2---(lh I < '^ith ... is equivalent to a system 

'Si,l7gi,l + Si,27gi,2 + 
^ S2,l7g2,l + S2,27g2,2 + 

.Si,l7gi,l + Si,27gi,2 + ••• + 0'qi,nnqi,nl = 

The questions are: what is the geometry underlying this parametrization? What is known 
about these sort of manifolds? What other properties of the resonance manifold are defined by 
a given summation rule? What additional information about these manifolds gives us the fact 
that they have many (often infinitely many) integer points? 

Another group of questions concerns transcendental dispersion functions. All our examples 
were constructed for rational and irrational dispersion functions and the theoretical results were 
based on some classical theorems on the linear independence of some sets of algebraic numbers. 
In the case of a transcendental dispersion function like = tanhafc similar reasoning can 
be carried out using the theorem on the linear independence of the exponents but it is not 
done yet. The question about special functions in this context is completely unexplored though 
very important. For instance, a dispersion function for capillary waves in a circle domain is 
described by Bessel function. Any results on their resonant interactions will shed some light 
on the nature of Faraday instability. 

One of the most interesting questions about the resonance manifolds would be to study their 
invariants, i.e. some new function / = /(m, n, 0) such that 0i + 02 = 03 implies /i + /2 = /s- 
An example of this sort of analysis is given in [15] for 3-wave interactions of drift waves with 
= ax/{l + y"^) but for real-valued wave vectors, x,y G R. Existence of the invariants is 
important because it is directly connected with the integrability of corresponding nonlinear 
PDE. Coming back to the physical language this means that the wave system possesses some 
additional conservation law. 

A very important task would be to develop fast algorithms to compute integer points on 
the resonance manifolds. The parametrization property allows to construct specific algorithms 
for a given dispersion function as it was done in [12] for 4-wave interactions of gravity waves, 
= (m^ + n?Y^^ . The work on the generic algorithm for a dispersion function = 0(k) 
( k = \/ TV? + or k = a/ + n"^ + P with integer m, n, I ) is on the way [T6] but it does 
not cover even simple cases like = m/n^ not mentioning transcendental dispersion functions. 
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